Quantum irreversibility, perturbation independent decay, and the parametric theory 

of the local density of states 
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The idea of perturbation independent decay (PID) has appeared in the context of survival- 
probability studies, and lately has emerged in the context of quantum irreversibility studies. In 
both cases the PID reflects the Lyapunov instability of the underlying semiclassical dynamics, and 
it can be distinguished from the Wigner-type decay that holds in the perturbative regime. The 
theory of the survival probability is manifestly related to the parametric theory of the local density 
of states (LDOS). In contrast to that the physics of quantum irreversibility requires to take into 
account subtle cross correlations which are not captured by the LDOS alone. 



I. INTRODUCTION 

The study of quantum irreversibility llj has become 
lately of much interest ]2|, ||, ||, ||, [| 0T|| due to its 
potential relevance to quantum computing, and to the 
theory of dephasing || [To| . Following JlJ we define in 
Section 2 the main object of the present Paper: This is 
the " fidelity" , also known as " Loschmidt echo" , which 
constitutes a measure for quantum irreversibility. 

The analysis of "fidelity" necessitates a generalization 
of the theory regarding "survival probability" pj}| . In 
Section 3 we remind the reader that the latter reduces 
to the analysis of the local density of states (LDOS) Jl2"| . 
Is it possible to make a similar reduction in case of the 
"fidelity"? At first sight such reduction looks feasible 
because the general physical picture looks very similar 
(Sections 5-6). 

In the present Paper we claim (Section 6), and prove 
by a numerical example (Sections 7-8), that the study of 
fidelity cannot be reduced to analysis of LDOS functions. 
Rather, it is essential to take into account subtle cross 
correlations which are not captured by the LDOS alone. 

The object of the present study is common to almost 
all the "quantum chaos" studies. Namely, to figure out 
what is the role of semiclassical mechanics in quantum 
mechanics. Whenever we find such (semiclassical) "fin- 
gerprints", we call them " non- universal effects". Most 
of the studies in the " quantum chaos" literature during 
the last 20 years have been devoted to figuring out the 
non-universal features of the energy spectrum. The main 
tool in singling out such features is a comparison with the 
predictions of random matrix theory (RMT). 

In the present Paper we use the same " philosophy" . 
Namely, we identify "non universal effects" by making a 
comparison with a corresponding random matrix model. 
On the other hand we discuss (Sections 10-11) a unify- 
ing theoretical picture that put the study of "quantum 
irreversibility" in the larger context of phase-space based 
semiclassical approach. 

An important ingredient in the understanding of " non 
universal features" follows from studies of the the clash 
between "perturbation theory", semiclassical theory, and 



RMT |h], Ell . A major realization is that semiclas- 
sical theory and RMT lead to different non-perturbative 
limits. Hence the resolution of the clash between the 
different theories involves the identification of different 
regimes of behavior. This applies i n g eneral to the anal- 
ysis of time dependent dynamics fll3f , and in particular 
to the analysis of wavepacket dynamics, decay of the sur- 
vival probability, structure of the LDOS Q, and natu- 
rally also to "quantum irreversibility" studies. 

Specifically, we distinguish in the present Paper be- 
tween regimes of "perturbative" and " non- universal" be- 
havior, and we define and study a billiard related model, 
where we have full control over the "borders" between 
these regimes. The conclusions are summarized in Sec- 
tion 12. 



II. THE "FIDELITY" 

Consider a system whose evolution is governed by the 
chaotic Hamiltonian 



H = H(Q 1 P;x) 



(1) 



where (Q,P) is a set of canonical coordinates, and x is 
a parameter. Later (Section 7) we are going to consider, 
as an example, a billiard system, where (Q, P) are the 
position and the momentum of a particle, while x is used 
in order to parametrize the shape of the billiard. Specif- 
ically, for a stadium we define x as the length of the 
straight edge, and adjust the radius parameter such that 
the total area is kept constant. 

Consider some Ho — TC(Q, P;xq), and define Sx = 
x — xa- Assume that 8x is classically small, so that both 
7^0 and TL generate classically chaotic dynamics of simi- 
lar nature. Physically, going from TLo to TL may signify 
a small change of an external field. In case of billiard 
system Sx parametrizes the displacement of the walls. 
Given a preparation 'J'o, the fidelity is defined as p3: 



M(t; Sx) 
m(t; Sx) 



\m(t; Sx)\ 2 

(* | ex.p(+mt) exp(~m t)\^o) 



(2) 
(3) 
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If "to is an eigenstate \E ) of Ho, then M(t; Sx) is equal 
to the survival probability P(t; Sx), which is defined as 

P(t;Sx) = \c(t;Sx)\ 2 (4) 
c{t;6x) = {E \exp{-int)\E ) (5) 

In the general case the preparation "So does not have to 
be an eigenstate of Ho- To be specific one assumes that 
is a Gaussian wavepacket. It is now possible to define 
a different type of survival probability as follows: 

P(i;wpk) = |c(i;wpk)| 2 (6) 
c(i;wpk) = (tfo|eJcpH«t)|tfo) (?) 

We assume that Sx is small enough so that we do not have 
to distinguish between H and Ho in the latter definition. 

One may regard Wo as an eigenstate of some prepara- 
tion Hamiltonian 7i wpk . Specifically, if Wo is a Gaussian 
wavepacket, then it is the groundstate of Hamiltonian of 
the type (P — Po) 2 + (Q — Qo) 2 that differs enormously 
from H. Thus we have in the general case three Hamil- 
tonians: 

• The preparation Hamiltonian W wpk 

• The unperturbed evolution Hamiltonian Ho 

• The perturbed evolution Hamiltonian H 

Above we have distinguished between two cases: The 
relatively simple case where H wpk = Ho- And the 
more general case, where we assume that the difference 
\\H„ pk — Wo 1 1 is in fact much larger compared with the 
perturbation \\H— Hq\\. The strength of the perturbation 
is controlled by the parameter 6x. 

III. THE LDOS FUNCTIONS 

Consider first the special case where 7"C pk — Ho- 
In such case the fidelity amplitude m(t; 5x) is just the 
Fourier transform of the local density of states (LDOS): 

p(u;5x) = Y,\Hx)\Eo)\ 2 S(iu - (Ejx) - E )) (8) 

n 

For technical reasons, we would like to assume that there 
is an implicit average over the reference state \Eq). This 
will enable a meaningful comparison with the more gen- 
eral case which is discussed below. 

In the general case, where H„ pk ^ Ho, one should rec- 
ognize the need in defining an additional LDOS function: 

p(^;wpk) =^|(n|* }| 2 8{lo - (E n - Eo)) (9) 

n 

In this context Eo is consistently re-defined as the mean 
energy of the wavepacket. Recall again that Sx is as- 
sumed to be small enough, so that we do not have to 
distinguish between H and Hq in the latter definition. 

The Fourier transform of p(aj;wpk) equals (up to a 
phase factor) to the survival amplitude c(t; wpk) of the 



wavepacket. The physics of c(t; wpk) is assumed to be of 
" semiclassical" type. We shall define what do we mean 
by "semiclassical" later on. The same notion is going to 
be used regarding c(t; Sx) if the perturbation (Sx) is large 
enough. 



IV. DEFINITIONS OF F AND 7 

In this Paper we measure the "width" of the LDOS of 
Eq.(3) in energy units Q, and denote it by T(Sx). A 
practical numerical definition of T(Sx) is as the width of 
the central region that contains 70% of the probability. 
This corresponds to the notion of " core width" in [jl2| . If 
Sx is not too large (see definition of " Wigner regime" in 
the next section), one observes that 

T(5x) oc Sx 

2/(1+9) (10 ) 

The value g ~ applies for strong chaos |L2| |l6|] , and it 
is the same as in Wigner's random matrix theory (RMT) 
p5| . In general (eg see Appendix) we can have Q < g < 1. 
In fact, the value g ~ 1 applies to our numerical model, 
which will be defined in Section 7. 

The decay rate of either the fidelity or of the survival 
probability (depending on the context) is denoted by 
'y(Sx). The semiclassical value of the decay rate, which 
is determined via a "wavepacket dynamics" phase-space 
picture Jll| , is denoted by 7 sol . The Lyapunov exponent 
is denoted by % t . 

In order to determine 'y(Sx) numerically one should 
plot M(t; Sx) against t, for a range of Sx values. In Sec- 
tion 7 we are going to define some model Hamiltonians 
for which we have done simulations. These are called 

• LBH: Linearized Billiard Hamiltonian 

• RLBH: Randomized version of LBH. 

• MBH: Modified Billiard Hamiltonian 

• RMBH: Randomized version of MBH. 

Fig. la displays the results of the MBH simulations. We 
see that the MBH decay is well approximated by expo- 
nential function (Fig. la). The dependence of the decay 
rate 7 M bh on Sx is presented in Fig. 2. The RMBH decay 
(Fig. Ib) is badly approximated by exponential function, 
but in order to make a comparison we still fit it to expo- 
nential. This is done in order to have quantitative mea- 
sure for the decay time. Thus we have also 7 RM bh(<5:e)- 

In Fig. 2 we also plot the LDOS width T(Sx) as a func- 
tion of Sx for the two models. As far as T(Sx) is con- 
cerned the two models are practically indistinguishable. 
The inset contains plots of T(Sx) for the other two mod- 
els (LBH, RLBH). In later sections we shall discuss the 
significance of the observed numerical results. 
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V. THE DECAY OF P(t; 5x) 



VI. THE DECAY OF M(t; Sx) 



The theory of the survival amplitude is on firm grounds 
thanks to the fact that it is the Fourier transform of 
the LDOS. According to |L2| there are three generic Sx 
regimes of behavior: 

• The standard perturbative regime. 

• The Wigner (or Fermi Golden Rule) regime. 

• The non-universal (semiclassical) regime. 

In the standard perturbative regime (Sx <C 8x c ) the 
LDOS function Eq.(|8[) is predominantly a Kronecker 
delta. This characterization constitutes a definition of 
this regime. For estimate of 5x c in case of billiards 
see Appendix. The survival amplitude is obtained via 
a Fourier transform of the Kronecker delta dominated 
LDOS function. This leads to a non-averaged m(t; Sx) 
that does not decay. On the other hand the Eq aver- 
aged m (t; Sx) has a Gaussian decay. The latter follows 
from the observation |Q that the first order correction 
Eq(x) — Eq(xq) has typically a Gaussian distribution. 

For intermediate values of Sx the decay of P(t; Sx) is 
typically of exponential type with 



7 = T(5x)/Ti 



(11) 



This is known as Wigner-type (or as Fermi golden rule) 
decay. It is a reflection of the Lorentzian-like line shape 
of the LDOS function. However, for large Sx we get into a 
non-universal (semiclassical) regime, where we can apply 
the " wavepacket dynamics" picture of (jll| . Thus we find 
a semiclassical decay with 



7 = 7.. 



(12) 



The Wigner regime, where Eq. (|Tl]) holds, is determined 
p2[ by the condition 



Y(8x) < h%, 



(13) 



This inequality can be re- written as Sx < Sx NV . The 
elimination defines a non-universal (system specific) 
parametric scale Sx NU . 

In the non-universal regime the width of the LDOS 
is semiclassically determined G|. In typical cases the 
width of the LDOS is proportional to the strength of the 
perturbation, hence 



7 BCl cx Sx 



(14) 



But in some exceptional cases % cl becomes perturba- 
tion independent. Specifically, for billiard systems l/7 sol 
roughly equals to the mean time between collisions, so 
we can write 



7s, 



7c: 



(15) 



It is important to realize that the perturbation indepen- 
dent decay (PLD) of c(t; Sx) in billiard systems, is a reflec- 
tion of the Sx independence of LDOS function p(u>; Sx) in 
the non-universal regime. See [16[ for a numerical study. 



A mature theory of fidelity is still lacking. However, it 
has been realized in ^ || that the same physical picture 
as in |Q arises: For very small Sx we have Gaussian 
decay [which corresponds to the Eq averaged decay of 
the survival amplitude]. For intermediate values of Sx 
we have Wigner-type decay with 7 = T(Sx)/H. For large 
Sx we enter into the semiclassical regime where one finds 
"Lyapunov decay" (25) with 7 w 7 scl w 7 cl . 

In complete analogy with the case of survival proba- 
bility studies we can define (via Eq.flL3)) an analogous 
parametric scale || that will be denoted by (Sa; NUD . The 
semiclassical value (7 Bd ) of 7 is not necessarily the same 
for P(t;Sx) and for M(t;Sx). Therefore in general Sx NU 
and <5x NUD are not necessarily identical. 

For simple shaped billiard system the 7 sc j of the sur- 
vival probability, the j sc i of the fidelity, and the Lya- 
punov exponent 7^ are all equal to the inverse of the 
mean collision time. The perturbation parameter Sx is 
defined as the displacement of the billiard wall. In the 
Appendix we derive the following result: 



(16) 



where 2ir/k is the De-Broglie wavelength of a particle 
with mass to, corresponding to the kinetic energy Eq = 
(Tik) 2 /2m. This results holds for hard walled billiard. 

In the following we want to demonstrate the distinc- 
tion between <5a: NUD and Sx NV . Therefore we consider a 
modified billiard Hamiltonian (MBH) for which 



^NUD Sx^y 



(17) 



The above inequality reflects the general case, in which 
the jsd of M (t; Sx) is different (smaller) from the 7 SC ; of 
P(t;5x). 

The fact that P(t;Sx) is a special case of M(t;Sx), 
and the fact that similar ideas (semiclassical decay ver- 
sus Wigner-type decay) have emerged in the latter case, 
naturally suggests that the same physics is concerned. If 
it were really the "same physics", it would imply that 
the main features of M(t; Sx) are determined by a simple 
minded theory that involves the LDOS function p(uj; Sx) 
in some combination with the LDOS function p(oj; wpk). 
It is the purpose of the following sections to demonstrate 
that a simple-minded theory is not enough. The semiclas- 
sical PID in the case of M(t; Sx) necessitates a non-trivial 
extension of the LDOS parametric theory. 



VII. DEFINITION OF THE MODEL 

Our model Hamiltonian is the linearized billiard 
Hamiltonian (LBH) of a stadium system Q. It can be 
written as 



U = E + SxB 



(18) 
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Here E is the ordered diagonal matrix {E n (xo)}. The 
eigen-energies of the quarter stadium billiard, with 
straight edge xq = 1, have been determined numerically. 
The perturbation due to 5x deformation, is represented 
by the matrix B. Also this matrix has been determined 
numerically as explained in 0. 

In the following numerical study we have considered 
not the LBH, but rather a modified Billiard Hamiltonian 
(MBH), which is obtained from the LBH by the replace- 
ment 

B„ m h-v G(n - m) x B nm (19) 

where G{n — m) is a Gaussian cutoff function. This cor- 
responds physically to having soft walls (for explanation 
of this point see Appendix J of p7t ). It is important to 
realize that the "exact" physical interpretation of either 
the LBH (as an approximation for the Billiard Hamilto- 
nian), or the MBH (as a soft wall version of LBH), is of 
no importance for the following. The LBH and the MBH 
are both mathematically "legitimate" Hamiltonians. 

In the next section we explain the numerical strategy 
which we use in order to prove our main point. This 
incorporates the random matrix theory (RMT) strategy 
which has been applied in [H in order to demonstrate 
that the semiclassical theory and RMT lead to different 
non-perturbative limits. The randomized LBH (RLBH) 
is obtained by sign-randomization of the off-diagonal el- 
ements of the B matrix: 

B m „ i-> ±B„ m (random sign) (20) 

The randomized MBH (RMBH) is similarly defined. 
The purpose in making comparison with "randomized" 
Hamiltonian, is the ability to distinguish between "uni- 
versal" and " non- universal" effects. Making such distinc- 
tion is a central theme in the "quantum chaos" literature. 
Usually such "comparisons" are made in the context of 
spectral statistics analysis, while here, following Q we 
are doing this "comparison" in the context of quantum 
dynamics analysis. 

VIII. THE NUMERICAL STUDY 

The first step of the numerics is to calculate the 
"width" T{5x), of the LDOS function p(ui;Sx). We know 
from previous studies jlj, |l6) that for hard-walled bil- 
liard system T(5x) shows semiclassical saturation for 
Sx > <5x NU , where 5 x NlJ roughly equals to De-Broglie 
wavelength (Eq.(|l^)). This implies PID for the survival 
probability. With the LBH we still see (inset of Fig. 2) 
a reminisces of this saturation. Note that k ~ 50 and 
hence Sx NU ~ 0.1. In contrast to that, with the RLBH 
there is no indication for saturation. This implies that 
non-trivial correlations of off-diagonal elements play an 
essential role in the parametric evolution of the LDOS. 
(See |l^] regarding terminology). 

By modifying the billiard Hamiltonian we are able to 
construct an artificial model Hamiltonian (MBH) where 



the two parametric scales are well separated (Sxyim <C 
Sx NU ). Thus within a large intermediate 8x range |20[| we 
do not have PID for P(t;Sx), but we still find PID for 
M(t;Sx). See Fig.2. 

In order to prove that the observed PID is not a triv- 
ial reflection of p(uj] wpk) we have defined the associated 
"randomized" Hamiltonian (RMBH). The LDOS func- 
tions (g) and d|) are practically not affected by the sign- 
randomization procedure: the sign-randomization proce- 
dure has almost no effect on T(Sx). In spite of this fact 
we find that the previously observed PID of M(t; 6x) goes 
away: we see (Fig.2) that for the MBH there is no longer 
PID in the relevant Sx range |^0| . This indicates that the 
PID was of semiclassical "off-diagonal" origin. 

We see that both qualitatively and quantitatively 
the sign-randomization procedure has a big effect on 
M(t;Sx). Therefore, we must conclude that the corre- 
lations of the off-diagonal terms is still important for the 
physics of M(t; 8x). This holds in spite of the fact that 
the same off-diagonal correlations are not important for 
the LDOS structure. This implies that the theory of 
M(t; Sx) necessitates a non-trivial extension of the para- 
metric LDOS theory. 

IX. THE SIMPLE MINDED THEORY 

The purpose of the present section is to explain what 
type of "fidelity physics" can be obtained if we do not 
take non-universal (semiclassical) features of the dynam- 
ics into account. Such theory is expected to be valid 
in case of RMT models. Let p et t(u); Sx) be the Fourier 
transform of m(t; Sx). It can be written as 

PeH {u>;5x) = Y,f{u')Ku-u') (21) 

where the summation is over energy differences uj' = 
(E n (x) — E m (xo)), and f(uj') is a product of the overlaps 
{n(x)\m(xo)) , and {m(xo)\^o) an d (^o\n(x)). It is clear 
that f(ui') satisfies the sum rule f(uj) — 1. On the 
other hand, if the number of principle components (par- 
ticipation ratio) of the LDOS is N, then the sum over 
|/(w)| gives N 1 / 2 . Thus we conclude that f(uj) should 
have random-like phase (or random-like sign) character. 
Therefore, if we ignore the system specific features, we 
can regard f(uj) as the Fourier components of a noisy 
signal. These Fourier components satisfy 

</M> = (22) 
(|/H| 2 ) = p(uj;wpk) xp(uj;Sx) (23) 

where p is, up to normalization, the auto convolution 
of p(w;wpk), and therefore equals to the Fourier trans- 
form of P(i;wpk), and has roughly the same width as 
p(uj] wpk). 

It is worth noticing that for Lorentzian line shape, 
which in general is not necessarily the case, Eq.(|2^) im- 
plies that m(t) is characterized by exponential correla- 
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tions with decay constant Y/2. This leads to decay con- 
stant r for M(t). The deviation of j(Sx) form T(Sx) in 
the MBH case can not be explained by Eq.(p3|), since the 
latter does not distinguish between the MBH model and 
the associated RMBH model. In order to explain the PID 
in the MBH case it is essential to take into account the 
non-universal (semiclassical) features of the dynamics. 



X. THE SEMICLASSICAL THEORY FOR P(t) 

The semiclassical theory of the survival probability is 
described within the framework of wavepacket dynamics 
in Ref . O] . The short time decay of c(£;wpk) reflects 
the loss of overlap between the initial and the evolv- 
ing wavepackets. On the other hand, due to the (in- 
evitable) proximity to periodic orbits, the survival am- 
plitude c(£;wpk) have recurrences. However, because of 
the (transverse) instability of the classical motion these 
recurrences are not complete. Consequently the long- 
time decay may be characterized by the Lyapunov expo- 
nent 7 C ; . Possibly, this " Lyapunov decay" is the simplest 
example for PID. It is PID because the size of the per- 
turbation [\\7i — 7iwpk||) is not relevant here. 

The semiclassical behavior of the survival probability 
has reflection in the LDOS structure. A relatively slow 
" Lyapunov decay" (due to recurrences) implies that the 
LDOS is "scared" @. Thus the semiclassical LDOS has 
a 'landscape' which is characterized by the energy scale 
hjd. Note that "scarring", in the mesoscopic physics 
terminology is called "weak localization" effect. 

The above semiclassical picture regarding c(i; wpk) can 
be extended |L2|, [l^, to the case of c(t;5x), pro- 
vided Sx is large enough. In the other limit, where 5x 
is small, we should be able to use perturbation theory 
in order to predict the decay rate. Thus we have here 
a clash of two possibilities: Having Wigner type decay 
with 7 = T(Sx)/h, or having non-universal decay (NUD) 
that reflects the semiclassical wavepacket dynamics. 

The crossover from the perturbative to the semiclas- 
sical regime can be analyzed jl2|, by looking on 
the parametric evolution of p{u>;5x). Depending on 5x 
the LDOS p(iv;Sx) has (in order of increasing perturba- 
tion) either standard perturbative structure, or core-tail 
(Lorentzian-like) structure, or purely non-perturbative 
structure [^]. The width r(<5x) of the "core" defines 
a 'window' through which we can view the semiclassical 
'landscape'. This landscape is typically characterized by 
h~"fd features, where 7 C ; is related to the classical dynam- 
ics. As Sx becomes larger, this 'window' becomes wider, 
and eventually some of the semiclassical landscape is ex- 
posed. Then we say that the LDOS contains a "non- 
universal" component [B2j. 



XI. THE SEMICLASSICAL THEORY FOR M(t) 

Whereas Lyapunov decay for c(t; Sx) is typically a 
"weak" feature [this is true for generic systems, whereas 
billiard systems constitute an exception], it is not so for 
m(t;Sx). By definition the trajectory of the wavepacket 
is reversed, and therefore the short-time decay due to a 
loss of wavepacket overlap is avoided. As a results the 
perturbation independent "Lyapunov decay" becomes a 
predominant feature (that does not depend on "recur- 
rences" ) . This Lyapunov PID has been discussed in Q . 

It is clear however that for small Sx we can use per- 
turbation theory in order to predict the decay rate of 
M(t;Sx). The question that naturally arise, in complete 
analogy to the P(t\ 5x) case, is how to determine the bor- 
der Sx NVU between the perturbative regime (where we 
have Wigner type decay) and the semiclassical regime 
(where we have NUD). 

The natural identification of Sx NVU is as the Sx for 
which T(Sx) becomes equal to h^sd. How 7 SC ; is deter- 
mined? There are two " mechanisms" that are responsible 
to the loss of wavepacket overlap. One is indeed related 
to the instability of the classical motion, while the other 
is related to the overall energy width of the wavepacket. 

The survival probability P(t; 5x) can be regarded as a 
special case of M(t;Sx), where the overall energy width 
of the wavepacket is the predominant limiting factor in 
the decay, the separation between the energy surfaces 
of TL and of TL$ is proportional to Sx. Consequently we 
typically have 7 SC ; oc Sx. 

In the prevailing studies of M(t; Sx), one assumes wide 
Gaussian wavepackets. Therefore the separation between 
the energy surfaces does not play a major role in the 
semiclassical analysis. Rather it is the instability of the 
classical motion that is the predominant limiting factor 
in the decay. Therefore one typically expects to have 
Isci ~ 7c;, which is independent of Sx. 



XII. CONCLUSIONS 

The above discussed criterion for the identification of 
the non-universal regime is in the spirit of spectral statis- 
tics studies p3fl . In the latter context it is well known 
that RMT considerations dominate the sub fi^d energy 
scale, while non-universal corrections dominate the larger 
energy scales. 

In the present Paper we have identified the "non- 
universal" regime for a billiard related model (MBH). 
The border between the perturbative regime and the non- 
universal regime in the context of P(t; Sx) is fe NU , while 
in the context of M(t; Sx) it is fe NUD . 

The parametric scales Sx m and 5a; NUD are similarly 
defined, but there is an important distinction between 
them. The first parametric scale marks the exposure of 
the semiclassical "landscape": either that of Eq.(||) or 
that of Eq.(EJ). The second parametric scale, as proved 
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by our numerical strategy, marks the exposure of cross 
correlations between the corresponding wave amplitudes. 



APPENDIX A: THE 5x NV FOR BILLIARDS 

The non-universal regime for billiard systems has been 
Here we would like to complete 
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missing steps in the generalization of this result. We use 
the same notations as in [|l2| [l6| . 

In the general case jl6| the bandprofile of the B matrix 
is determined by the semiclassical formula p4| 



IB, 



A 

2nh 



C 



E n — Er. 



(Al) 



where A oc l/k d ~ 2 is the mean level spacing, k is the 
wavenumber, and d — 2 is the dimensionality of the Bil- 
liard system. The power spectrum of the motion 



C(w) = const x k 3+9 /oj 9 



(A2) 



is the Fourier transform of a classical correlation function. 
Here g — corresponds to strong chaos assumptions, 
while < g < 1 is more appropriate for our type of 
system due to the bouncing ball effect. The width of the 
LDOS is determined using a procedure which is explained 
in (l^l, leading to Eq.(9) there. Namely, 



T(Sx) = A x {5x/Sx c ) 2/(1+9) 



(A3) 



where Sx c oc k 
Eq.(8) of @. 



-((l-ff)+(l+<7)d)/2 



is the generalization of 



Form Eq.(A3) it is clear that Sx c should be inter- 
preted as the deformation which is needed in order to mix 
neighboring levels. In the standard perturbative regime 
(Sx <C Sx c ) first order perturbation theory is valid as a 
global approximation. Otherwise, if Sx > Sx c , we should 
distinguish between a non-perturbative "core" of width 
r, and perturbative "tails" that lay outside of it. 



The expression for T can be re-written as 



T(Sx) » Ti lci x (kSx) 2 ^ 1 ^ 



(A4) 



where j c i oc k is roughly the inverse of the ballistic time. 
In our numerical analysis (7| we find that T sa 0.36fc 2 x Sx, 
corresponding to g = 1. The non universal scale Sx NU , 
as well as Sx NVU , are determined by the requirement 
T(5x) = h^d. Hence we get Eq.(|i~6|), which holds irre- 
spective of the g value. The latter claim has been stated 
in 1 12 without a proof. 
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FIG1: (a) The decay of M(t; Sx) in the MBH case, (b) The same 
for the randomized MBH (RMBH). We use dimcnsionless units of 
time that correspond to stadium billiard with straight edge xo = 1, 
particle with mass m = 1/2, wavenumber k ~ 50, and K = 1. 
The values of the perturbation strength are (from the top curve to 
bottom): Sx = 0.0125 * i with i = !,-■■, 11. 
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FIG2: The LDOS width T, and the decay constant 7 from the 
MBH/RMBH simulations. The dotted line is the classical Lya- 
punov exponent. The inset is T in the LBH/RLBH case. 



